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Abstract 

Coherent stochastic resonance is explained in terms of first passage time den- 
sity functions. Scahng relation between first passage time density functions at 
resonance for different lengths of the medium is obtained. A formula for first 
passage time density function at resonance is derived in terms of two univer- 
sal functions, which demonstrates the universal feature of coherent stochastic 
resonance. A simple approximate expression for first passage time density 
function at resonance is proposed which is shown to explain the behavior at 
resonance fairly accurately. 
KEY WORDS: Synchronization, coherent stochastic resonance, first passage time density 
functions 

PACS number(s):05.40.+j 



1 



Typeset using REVTgX 



I. INTRODUCTION 



There has been a large deal of interest in the understanding of mechanism of interplay 
between random noise and a deterministic periodic signal. It has been found that they act 
cooperatively in some nonlinear systems Nature, presumably, exploits this fact to its 

advantage to tune in to a desired signal [^0. Even for linear, spatially extended stochastic 
system the cooperative response between noise and deterministic periodic force has been 
demonstrated very recently 

The first passage time is a useful tool to investigate the diffusive transport property in 
a medium. The theory of first passage time has been worked out in great detail for both 
infinite medium and explicitly time- independent diffusive processes [PHTTH. However, for 
explicitly time-dependent processes and in finite medium an analytic closed form expressions 
are not available. In this respect also this problem attracts much attention to the scientific 
communities. 

The process that we are concerned is an overdamped system, having no intrinsic fre- 
quency, driven by a periodic force and embedded in a noisy environment which is taken to 
be Gaussian white for simplicity. The motion is constrained between two traps. The first 
analysis of this process has been done for a random walk on a lattice numerically and in 
a continuous medium with a periodic signal of small amplitude perturbatively [Q. Their 
results indicate that the oscillating field can create a form of coherent motion capable of 
reducing the first passage time by a significant amount. In order to investigate the reason for 
this cooperative behavior of random noise and deterministic periodic signal, this problem has 
been formulated in much simpler terms by approximating the sinusoidal periodic signal by 
the telegraph signal and later it has been concluded @] that the telegraph signal cannot pro- 
duce any cooperative behavior. Subsequently, it has been shown that if we approximate 
the sinusoidal signal by a multi-step periodic signal [explained below] one would recover the 
coherent motion. These authors |^ have also demonstrated why single-step telegraph signal 
cannot produce any coherent motion as conjectured by the earlier worker H] and multi-step 
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telegraph signal does. 

While some authors [0 demonstrated the non-monotonous behavior of the mobility as a 
function of forcing frequency by calculating the net mass leaving at one of the boundaries 
for weak force, others have demonstrated the non- monotonous behavior of the mean first 
passage time (MFPT) for small amplitude |^ and for arbitrary amplitude ||^ in a continuous 
medium. These authors have also demonstrated that the variance of the first passage 
time density function (FPTDF) also attains a minimum at the same frequency for which 
the MFPT shows a minimum when they are plotted as a function of forcing frequency. 
This happened for any length of the medium. Such an enhancement of the mobility at 
some definite forcing frequency depending on the length of the medium clearly exhibits the 
maximum cooperation between noise and periodic signal and is known in the literature as 
coherent stochastic resonance (CSR). 

All the previous analyses 0,^,0 concentrate mainly on the first order moment (MFPT) 
of the first passage time density function (FPTDF) except in p], where it has been shown 
that the FPTDF at resonance at different lengths of the medium overlap over each other 
when they are plotted as a function of the cycle of the periodic signal. This clearly indicates 
that the resonance feature is, in some sense, universal in nature. To investigate this feature 
further, we in this paper concentrate on the full profile of FPTDF and study the effect of 
synchronization between noise and deterministic periodic signal in terms of FPTDF as we 
change the frequency. It is shown that a single dominant peak of FPTDF occurs at resonance 
where the synchronization is maximum. This feature is quite distinct and is not observed 
when the system is away from the resonance. Further, a scaling relation is obtained between 
FPTDFs at resonance for different lengths so that one could map FPTDF at arbitrary length 
on a single FPTDF at our choice by such scaling relation. This feature together with the 
cycle variable of the periodic signal as argument of FPTDF conclude the universal feature 
of FPTDF at resonance. 

The paper is organised as follows. In Sec. II, we give the formulation of the problem. 
The results of the calculations based on the derived formulae are presented in Sec. III. Apart 
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from the explanation of synchronization in terms of FPTDF, the scahng behavior of FPTDF 
is also presented in this section. The FPTDF at resonance is expressed in terms of two 
universal functions, which clearly isolates its dependence on the length of the medium. 
Further, a simple expression of FPTDF at resonance is also obtained in the same section 
which is shown to be fairly accurate in explaining the behavior at resonance. Finally, a few 
concluding remarks have been added in Sec. IV. 



II. FORMULATION OF THE PROBLEM 

The system that we consider is described by a Fokker Planck equation 

= -As^nnt^-^ + D^\^f\ (1) 
ot ox ox^ 

where A and fl are the amplitude and frequency of the sinusoidal signal, and D is 

the strength of the Gaussian white noise. The motion is confined between two absorbing 

boundaries at a; = and x = L ; i.e., Eq.(|^) is supplemented with absorbing boundary 

conditions p{0, t) = p{L, t) = 0. In terms of the dimensionless variables 

C = iA/D)x, e = {A^/D)t, uj = Q/{AyD), (2) 

the Eq.(|l|) is conveniently written as 

The boundary conditions read in terms of the new variables (2) as p{0, 9) = p{A, 6) = 0, 
where A = {A/D)L. In the following we calculate all the physical quantities in terms of 
these new variables and if required, one may translate all the interpretations in terms of the 
usual variables by the transformation equations Eq.(0). 

No analytic solution exists for Eq.(^) with the boundary conditions as mentioned. We 
thus introduce a scheme to approximate the force sinu;6' as a multi-step periodic signal [§. 
This scheme is in contrast to the procedure adopted in 0, where they discretise the Eq.(^ 
using finite difference method and simulate the problem on a lattice of space and time. 



The construction of multi-step periodic signal is as follows. We divide the half cycle 
of the signal by {2p + 1) intervals so that each interval in the horizontal ^-axis is of size 
{A9/{2p + 1)) with uj A 6 = n . We define {2p + l)numbers Sk along the vertical axis as 

Sk = ^ ; ^ = 1, 2, (4a) 

Sp+i = 1 (4b) 
"5p+i+r = "Sp+i-r ; r = 1, 2, p. (4c) 

Each number si. is asso dated with the interval ^^jf^ <0 <^ with k = 1,2, {2p + 1). 
The Eq.(^ clearly shows that 

< Si < S2 < ... < Sp < Sp+i = 1 > Sp+2 > Sp+3 > ••• > ■S2P+1 > . (5) 

Eq.(^ states that in order to reach the maximum value (= 1) of the signal from the zero 
level we have to have {p + 1) step up and from the maximum to the zero level we have 
{p + 1) step down. This is for the positive half-cycle. For the negative half-cycle similar 
constructions have been done with the replacement — ^ —Sk,^k and each number —Sk is 
associated with the interval A0[1 + ^] < 6 < A0[1 + ^] with k = 1,2, {2p + 1). 
This aproximation for the full one cycle of the sinusoidal signal (as shown in Fig.l) is then 
repeated for the next successive cycles. 

One may however note that the uj which we have defined for this approximated signal 
is not the same as that of sinusoidal signal because the Fourier transform of sinusoidal 
signal would give only one frequency while this approximated signal in the Fourier space 
corresponds to many sinusoidal frequencies specially because of its sharp discontinuities. Yet 
we urge this approximation because in each interval the equation become time- independent. 

In the future development we associate the index n for the positive half-cycle and index m 
for the negative. Index i will refer the cycle number. Since the Fokker-Planck equation Eq. (^ 
in each interval will be that for a constant bias, we can express the conditional probability 
density function p{^,6 \ ^',0') in terms of complete orthonormal set of eigenfunctions m„(0 
satisfying the boundary conditions 'U„(0) = ^^(A) = 0. 



n 

where 

1 TlTTc 

n±(0 = (2/A)5ea;p(±se/2)szn-^ , (7a) 

^^^^ N 

An = ^ + J . (7b) 

with s as the corresponding value of Sk in the appropriate interval where the conditional 
probability is being decomposed. The conditional probability density function in any inter- 
val, say /, can then be calculated from the previous history by convoluting it in each previous 
intervals: 

Pi^i, Oi I 6, 0i)= f- f d^i-id^i-2...d^2 ri Pi^v (^j I ■ (8) 

i=2 

For the negative half-cycle the calculation of probability density function is similar except 
that we have to replace the index n by m and the probability density function is decomposed 

as 

I 6') = Y: u^^iOuiiOexpi-XUO - e')] , (9) 

m 

where the expressions for m^(0 are same as in Eqs.(^. 

The survival probability at time 6 when the particle is known to start from ^ = at 
^ = is defined as 

S{9\^o)= r«e,^Uo,0) . (10) 
Jo 

The first passage time density function (FPTDF) g{6) is defined as 

9(« I &) = . (11) 

Physically, g{6)d6 gives the probability that the particle arrives at any one of the boundaries 
in the time interval 6 and 6 + d9. From this density function one can calculate various 
moments: 
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< >= / dee^g{e) . (12) 





From Eg . ([T^) one can easily calculate mean first passage time(MFPT) < 6 > and the 
variance o"^ =< 6^ > — < 6 >^ of the density function g{6). 

It is then quite straightforward to calculate the FPTDF at any interval of any cycle. We 
will write down the final formulae: 

g+{9 I 6) = X A„^,^^^,(^_^,^^ X exp{-A„(,^^^,(^_^,^^ x - 2{z - 1) A 9]} 

X-Fj-i('U„^2p+i)(,_i)+i) 

■,2{t-l)A9<9< {2{i - 1) + A 9 , (13a) 

2p+l 

q+{9 I fn) = X Art„ i,^-,, X expi-A^,, - 2(i - 1) A ^11 x 

"{2p+l)(i-l)+0+l) I n(2p+i)(i_i)+j J 



A^ fc-i 



;(2(.-l) + ^)A.<.<(2(.-l) + |±l))A. 

;A; = l,2,...,(2p-l) , (13b) 



I eo) = X A„,,^^,,^ X exp{-A„,,^^,,J^ - (2. - 1) A 9]} 



X 



^ ('"n(2p+i)»' ^"(2p+i){i-i)+i'' ^ -^«-l'^'"n(2p+l){«-l)+i) 

(2(i - 1) + ^^^) A^ < ^ < (22 - 1) A^ , (13c) 



I ^o) = C„^2p+i)(,-i)+i ^ ^-,2p+i)(,-i)+i X exp{-\m^,^^^^,,_,,_AO - (2^ - 1) A 9]} 

™{2p+l)(i-l) + l I "(2p+l)i 
X^^(M„(2p+i)i5 '"n{2p+i)(i_i) + i) X -Fj-l('"n(2p+i)(i_i) + i) 

; (22 - 1) A < < ((22 - 1) + -^) A 9 , (13d) 

2p + 1 
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g40 I Co) = C'^(,,+,„,_i)+(,+,) X A^(,^^,,(,_,,^(,^,, X exp{-Xm^^^^^^^,_,^^^,^,^ [9 - {2i - 1) A 9]} x 

k 

m(2p+l)(i-l) + 0+l) I m(2p+l)(i_l)+3 J 

xexp[2^ + 1 ^^'^'"(2p+i)(i-i)+(fc+i) ~ X^''^'"(2p+i)(i-i)+o+i))] 

X <^ I ^ 

^ ^ "'m(2p+i)(i_i)+i I "'n(2p+i)i ^ 

X^^(^in(2p+l)i'^n(2p+l)(i_l)+l) X -^j-l('"n(2p+i)(i_i)+i) 

;* = l,2,...,(2p-l) , (13e) 



where 



9^0 I Co) = X A^^,^^,,, X exp[-X^^^^^^^^{e -2iA 9)] 

^ "»(2p+i)i' "^(2p+i)(i-i)+i^ "'m(2p+i)(j_i)+i I "'n(2p+i)j ''^ 

X^"^(lin(2p+i)i' ^•''n(2p+i)(i_i)+i) ^ -^«-l('"n(2p+i)(i_i)+i) 

; ((2^ - 1) + Ae<e<2iAe, (i3f) 



C+= f'^diuiii) , (14a) 
Jo 



Cm 



tdiu-{0, (14b) 
Jo 



^'^('"n(2p+i)i''"n(2p+i)(i_i)+i) (r, , ■i)-^n(2p+l){i-l)+l 



Ag 

m< li~ I m"*" > 

'l(2p+l){«-l) + (i + l) I n(2p+i)(i_i)+j 

J = l 



A0 

xe^^^[-(^^)\2.+i)(.-i)+o+i)]} > (14c) 



^ ('"m(2p+i)i'^m(2p+i)(i_i)+i) ^^^'I (2^.^ ]^)'^"^(2p+i)(«-i)+i] 



2p 

m< I > 

m(2p+l)(i-l) + 0-+l) I m(2p+l)(i_l)+3 

A^ 

X exp[- ( 2^ _^ ^ ) Ar„(2p+i)(i-i)+o+i) ]} (14d) 
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and the functions Fi are generated through the recursion relation: 

FAu~ ) = < u~ \ u~ > xA~(ut, ,u~ ) 

"(2p+l)i+l^ "(2p+l)i+l I »^(2p+l)i ^ "^(2p+l)i' "^(2p+l)(i-l) + l 

X < I > xA^(u^ ) 

™-(2p+l)(i-l) + l I "(2p+l)i ^ "(2^+1)1' »^(2p+l){i-l) + l' 

x^i-i(«« ) , (15) 

* -^^ "(2p+l)(i-l) + l^ ' 

with Fo{u~J = u~_^{$,o). The angular bracket in any equation implies dot product of the 
corresponding functions, for e.g., 

<u+\u->= [^d^u+iOu~iO ■ (16) 







The cycle variable i runs over positive integers;i.e.,z = 1,2,3, .... The positive and negative 
symbols of the FPTDF indicate its values over positive and negative part of the cycles re- 
spectively. In all these expressions, viz., Eqs.(13)-([T5|), any subscript either n or m or both 
wherever they appear more than once the summation over them are implied. The effect of 
history is explicit in the expressions for FPTDF, in particular, through recursive relation 
between cycle variable i through Eq.(|T3p and through the product symbols in Eq.(13b), 
Eq.(13e), Eq.(14c) and Eq.(14d). The Eqs.(13) show how the non Markov character is sys- 
tematically incorporated through the convolution Eq.(|^) of the probability density functions. 
Once the FPTDF g{6 \ ^0) is obtained from these formluae, the various moments could be 
obtained by employing Eq.([I^). Evaluation of FPTDF and other relevant quantities requires 
sum of infinite series which must be truncated in order to obtain a final result. Convergence 
of FPTDF is ensured by gradually increasing the number of terms (i.e., number of eigenval- 
ues) for the calculation. The process is truncated when MFPT does not change up to two 
decimal point of accuracy with the change of number of terms. 

III. RESULTS AND DISCUSSIONS 

We have already demonstrated that p = approximation or usual telegraph signal 
approximation of the sinusoidal signal does not produce any resonance. In all the following 
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calculations we restrict to p = 2 , although p = 1 approximation of the signal does produce 
the resonance. But the p = 2 signal approximates better than the p = 1 signal. Further we 
take ^0 = A/2 , i.e., the particle is assumed to start from the mid point of the medium. 

A. Synchronization between noise and deterministic signal 

The FPTDF g{6; u) for different frequencies are plotted as a function of 6 for a given 
length A = 20 in Fig. 2. For very low frequency (for e.g.,ti; = 0.005), the particle is acted on 
by a constant force s = Si (in the multi-step periodic approximation) almost all the time 
before it reaches the boundary. The FPTDF curve shows that we have only one maximum 
in the entire 6'-regime. 

As the frequency slowly increases, the curve shows that apart from only one profile 
appearing left and showing a maximum like before, a small peak at large 6 also shows up 
(for e.g., uj = 0.01). This is because the particle at this frequency, although most of the 
time sees a constant bias s = si, at later time it might encounter a increased bias s = S2 
for a short time before it reaches the boundary. As J g{6)d6 = 1, therefore the area under 
the major profile decreases and this decrement has been compensated by the small peak 
at the right end of the profile. In other words, the g{6) for this frequency consists of two 
curves superposed over each other. Further increase of frequency (for e.g.,^^ = 0.02) shows 
the similar behavior as with uj = 0.01 but here the first peak area decreases and second peak 
area increases more than that with uj = 0.01. 

What we have been doing is that we follow the g{6) profile as we increase the frequency 
and observe a systematic change in the constituent profiles. As the frequency increases 
more, other biases Sk with k > 2 start playing their roles in constructing the FPTDF profile. 
For e.g.,ti; = 0.03,0.04 , we have another little peak at extreme right starts developing. As 
frequency further increases, this last peak starts developing further, while the first peak 
which has been dominant at very low frequency starts slowly decreasing. 

Our previous study with this scheme |^ shows that the MFPT has a minimum for this 
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length (A = 20) at = 0.1. We have identified this frequency as resonant frequency because 
we have argued that the synchronization between the signal and noise is maximum at this 
frequency causing the enhancement of probability of reaching the boundary at a short time. 
The FPTDF profile at resonance evidently shows a dominant single peak at this frequency. 
More smoother and nice single dominant peak at resonance is observed for higher length 
[see for e.g., Fig. 7]. Slightly below this resonance frequency {uj = 0.07, 0.08, A = 20) the 
FPTDF profile contains two peaks, the maxima of them are close to each other and the 
peak which has been dominant at very low frequency diminishes so much that only very 
close observation of its nature starting from low frequency would help identify its existence. 
We further note that the constituent profiles of FPTDF adjust themselves in a very distinct 
way, as the frequency increases, in order to produce a dominant maximum at the resonance. 

In order to explain synchronization we must interpret the cooperation in terms of cor- 
relation between the components of FPTDF profile. For this purpose we start looking at 
the FPTDF profiles for the frequencies higher than the resonant frequency. We have argued 
that at resonance, the synchronization becomes at its peak. After this resonant frequency, 
another small and distinct peak (for e.g., u = 0.13) starts developing. The area of the 
dominant peak starts reducing. As frequency increases further beyond the resonance fre- 
quency, more number of peaks start coming up over a little background (for e.g., u = 0.2). 
After the resonance, the peaks are quite distinct and identified separately over the visible 
background. Following the dominant peak beyond resonance, it is natural to identify the 
area of the peak as a degree of synchronization. Our previous study with this scheme Q 
shows that the MFPT increases beyond the resonance. In these figures we find that the 
dominant peak whose area is maximum at resonance decreases as we go beyond the reso- 
nant frequency. In other words, the degree of synchronization is getting reduced beyond 
the resonance frequency. Clearly we can find the areas of the individual distinct peaks and 
the background . The profile of the background is obtained by fitting the minima of the 
FPTDF profile. Going backwards from high frequency side it is clearly visible that the 
background area diminishes as we approach the resonance ; Most of the area is taken up 
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in producing the major dominant peak at resonance. The question naturally would be how 
one can distinguish the components, namely the background and the peaks as we decrease 
the frequency below the resonant frequency. 

The clue to the answer of this question lies in the FPTDF profile at very low frequency. 
As we have already noticed that at very low frequency we have only one peak (for e.g., 
oj = 0.005). The area of this curve gradually decreases as the frequency increases and 
approaches its minimum at resonance. After the resonant frequency, it appears clearly as 
background on top of which the peaks are occuring. It is evident that this curve does not 
really play a role in making the dominant peak at resonance. The other two peaks which 
appear at larger frequencies (say, at u = 0.07, 0.08) start adjusting themselves constructively 
to produce a large dominant peak at resonance {u = 0.1). Once we identify the background 
in the low frequency regime by following the component corresponding to the curve at 
u = 0.005, the area of the rest of the profile of FPTDF at any frequency over the background 
could be easily calculated. As we have argued that the two peaks really overlap each other 
especially near the resonance and adjust themselves constructively, the area under these 
should be a measure of synchronization before the resonance. 

Having identified the background, so far as the resonance is concerned, we try to fit 
this component in each FPTDF profile. As mentioned before, for very low frequency the 
full FPTDF gives the profile of the background. We have also noted that at extremely 
low frequency, in our multi-step periodic signal scheme, the particle experiences a constant 
bias s = Si throught the entire time. But no simple analytical formula of FPTDF in such 
situation exists. However, for semi-infinite domain, if we have an absorbing boundary at 
^ = A and the particle starts from ^ = A/2 the exact analytic expression for FPTDF was 



evaluated by Schrodinger |12|. The expression reads as 



9s{0) = ^^J^exp{-{A/2 - 3,6)^46} . (17) 

We plot this expression and superpose over actual FPTDF curve at u = 0.01. The plot is 
shown in Fig. 3. It appears that although for low 6 this function gives excellent fit, for large 
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6 it overestimates the background. Therefore we modify this function in the following way 
to arrive at a better fit of the background. 

g{e)=gs{e) ;O<0<^^o 

= gs{0)exp{-(3{d-dof} ■,0>0o (18) 

The fitted curve with 6q = 45.14,/? = 5.66 x 10~^ is shown in the same figure [Fig. 3(b)]. 
As frequency changes and approaches resonance, Oq and /? change. The area under this 
curve, Pf,, is plotted in Fig. 4a and found decreasing before resonance. After resonance, the 
profile of the visible background is obtained by joining the minima of FPTDF profile at 
each frequency. When the minima are plotted as a function of 6, it fits with the same curve 
Eq.(p!8D with 6q = and Si = 0. Only j3 is found changing slowly and the area under 
this curve is found increasing as the frequency increases beyond resonance (Fig. 4a). This 
concludes the characterization of the background. 

The area of the peaks over the background, P, responsible for constructing resonance 
before the resonance frequency is also plotted in the same figure (Fig. 4a) and found increasing 
with frequency. At resonance, as expected this area denoting the degree of synchronization 
attains a peak. As we follow the area of the dominant peak beyond resonance, it is found 
decreasing with frequency showing the degree of synchronization is getting reduced as we 
go beyond the resonant frequency. This is also demonstrated in Fig.4a. 

Further, beyond the resonance, as mentioned before, apart from the dominant peak 
there are other smaller peaks appearing at the top of the background. The areas, Pi,P2, 
under first two sucessive peaks beyond the dominant peak are also plotted as a function of 
frequency in the Fig.4b. These curves also show non-monotonous behavior. This shows that 
the synchronization might also occur at higher frequency but as the maximum area covered 
by these peaks are much less, they denote synchronization of lesser degree. This completes 
the explanation of synchronization behavior between noise and deterministic periodic signal 
as we change the frequency. 
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B. Scaling of FPTDF at resonance 



We have already noted that the FPTDF exihibits a dominant peak at resonance. This 
feature is observed for any length. In our previous study 0] we have found a simple relation 
between the resonant frequency uj* and the length of the medium A. The expression reads 

as 

uj* = 2/A . (19) 

It is further interesting to observe that the normalised FPTDF g{9)/uj* when plotted as 
a function cycle variable [co'*(A)0], they overlap over each other. The curves for different 
lengths are plotted in Fig. 5. This figure clearly shows that at resonance cycle variable is the 
correct argument to describe the resonance feature. We may further note that such scaling 
of the argument of FPTDF would not be possible for any frequency other than the resonant 
frequencies because any frequency which is not the resonant frequency for one length may 
turn out to be the resonant frequency for some other length and the features of FPTDF are 
different for resonant and off-resonant frequencies as has been observed from Fig. 2. 

Having found the proper scaling of the argument of FPTDF, it is natural to enquire 
whether the FPTDF, / = g{9)/uj* can also be scaled properly, so that once f{x) with 
X = uj*9 is found for one length, the function f{x) can be obtained for any arbitrary length. 

In order to investigate this issue, first we observe that going from /(a;;Ai) to /(a;;A2), 
where Ai and A2 are two different lengths we have to multiply one function by different 
amount depending on the value of the argument. Thus if at all any scaling relation exists, 
this should be of the form 

f,K{x) = /A(x)/i"(^) (20) 

It is indeed found to be true with a{x) given in Fig.6. We note that a{x) is universal. 
From the relation (20) it is obvious that 

A-°(-)/A(a;) = [/iA]""^"VM(^)- (21) 
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The Eq.(0) shows that the function A'"^^^ f/^{x) is independent of A and it is again universaL 
We verify this result in our calculation of FPTDF for different lengths and call this function 
U{x). The lnU{x) is also plotted in Fig. 6. Thus we could express fA{x) at resonance in 
terms of two universal functions U{x) and a{x) as 

/a(x) = f/(x)A"(^) (22) 

The expression (22) clearly isolates the dependence of FPTDF on length of the medium A 
at resonance. Therefore, for any arbitrary length A the FPTDF at rsonance can be found 
from this simple expression. 

We next fitted the universal functions a{x) and lnU{x) with the simple functions as 

a{x) = —aaipix) — Ca , (23a) 
lnU{x) = auip{x) — cu , (23b) 

where aa,auyCa,cu are positive constants and ip^x) = (1 — ■^)^/x with x = Xm being the 
position of the maximum of f{x). With these approximate expressions (23) the normalised 
FPTDF /a(x) at resonance takes a very neat expression 

U{x) = C( , (24) 

with Ao = 10.312, C = 0.5148. The normahsation of FPTDF suggests that J fA{x)dx = 1. 
It is indeed found that this approximate form (24) would result normalisation constant, for 
the lengths we consider, almost independent of A and very close to unity. The approximate 
form f^ix) is plotted for a typical length A = 40 and compared with actual data points 
in Fig. 7. It is giving an excellent agreement. Further with this form, the MFPT is found 
for different A. It also yields linear behaviour in agreement with our previous study [§. 
This shows that not only we have found a scaling behaviour (20) which takes f{x) from one 
length to the other, it is also possible to obtain a simple expression (24) for /A(a;) for any 
arbitrary length. 
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IV. CONCLUDING REMARKS 



We consider a diffusive transport process perturbed by a periodic signal in continuous 
one dimensional medium having two absorbing boundaries. No perturbation of the signal 
amplitude is assumed in this formulation. The scheme we have adopted in this formulation 
is different from the finite difference method to solve the partial differential equation Eq.(^ 
as adopted in The cooperative behaviour between the deterministic periodic signal and 
random noise leading to coherent motion is demonstrated in terms of FPTDF profile. We 
show explicitly that synchronization between the signal and random noise is maximum at 
resonance where FPTDF exhibits a major dominant peak. 

An important characteristic that we observe is that the cycle number is the correct ar- 
gument to describe FPTDF at resonance. We also show that there exists a scaling relation 
between FPTDF at various lengths through some universal function a{x). The exact expres- 
sion for FPTDF at resonance is obtained in terms of two universal functions which clearly 
isolates its dependence on the length of the medium. This is an universal feature of coherent 
stochastic resonance. Further a nice simple approximate form for FPTDF at resonance is 
obtained which is shown to be fairly accurate in explaining the behavior at resonance. This 
form may be of use in obtaining quicker result in cases where more complex situation has 
been called for. 

Although we restrict our calculation with a three-step {p = 2) periodic signal, the for- 
mulation is quite general and applicable for any approximation with an arbitrary number 
of steps. This formulation can also be applied to any arbitrary continuous periodic signal. 
For {p = 1) periodic signal the background could be identified as before. The resonance 
occurs for the length A = 20 exactly at the same frequency {uj = 0.1). We obtain a major 
dominant peak at resonance. The only difference that we observe is that the component 
of FPTDF just before resonance which makes the peak at resonance constitutes only one 
peak instead of two as with p = 2. The reason for which is tailored as each half-cycle of the 
sinusoidal signal is approximated with three steps for p = 1 as against with five steps with 
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p — 2. However, after resonance frequency the behaviour is similar as with p — 2 signal. 

The peak positions of normahsed FPTDFs at resonance occur very near to a quarter of 
a cycle. From Fig. 5 we observe a slight deviation of the peak positions but we believe that if 
the sinusoidal signal is approximated by more than three-step periodic signal, the position 
of all the peaks will be the same. 
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FIGURES 

FIG. 1. Sinusoidal signal(dashed curve) and approximated three-step(p=2) periodic signal(solid 
curve) for the full one cycle as a function of 6. 

FIG. 2. (a) FPTDF g{6) as a function of 9 for A = 20 before and on resonance for frequencies co 
= 0.005, 0.01, 0.02, 0.03, 0.04, 0.07, 0.08 and 0.1 (resonant) (b) FPTDF g{e) as a function of 9 for 
A = 20 after resonance for frequencies u = 0.13, 0.15, 0.2 and 0.3 respectively, \p = 2,^o = A/2]. 

FIG. 3. FPTDF g{9) as a function of 9 for frequency oj = 0.01 and A = 20(filled circle) 

[p = 2,^0 = A-/2] (a) The function gs{9) in Eq.(17) as a function of 9 (dashed curve) (b) The fitted 
function g{9) in Eq.(18) as a function of 9 (solid curve). 

FIG. 4. (a) The area of the background Pf, (filled square) and the area of the dominant peak 
P (filled circle) as a function of frequency uj (b) The areas Pi,P2 (filled circle) as a function of 
frequency lo. 

FIG. 5. The dominant peaks of g{9)/aj* at resonant frequencies for different lengths 
(A = 30,35,40,45,50) as a function of u>*9. The lowermost curve is for A = 30, and upper 
curves correspond to higher lengths respectively, [p = 2,^o = ^f^]- 

FIG. 6. The universal functions a{x),lnU{x) as a function of x. 

FIG. 7. (a) /a (a;) as a function of x for A = 40 at resonance (filled circle) [p = 2,^o = A/2] (b) 
/a(x) according to the formula (24) as a function of x for A = 40(solid curve). 
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